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Abstract 

The sodium reordering in Na^CoC^ in the vicinity of room temperature is rationalized at high x 
in terms of phase transitions between square and striped phases. A striking hexagon-of-hexagons 
diffraction pattern observed for £=0.78 can be reproduced using coexisting square and striped 
phases that are related by simple shear deformations. All compositions exhibit a partial melting 
transition to a disordered stripe phase just below room temperature, which alters the topology of 
the electrical conduction pathways. 

PACS numbers: 61.05.cp, 61.05.fm, 61.50.Ks 
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Nanopatterning of sodium vacancy clusters in sodium cobaltate, Na^OoC^, templates 
the Coulomb landscape on the Co layers and, therefore, controls the electronic and mag- 
netic properties^. The concentration of sodium, x, can be altered electrochemically, directly 
changing the number of electrons in the triangular Co layers^. Furthermore, for a fixed 
composition the physical properties depend sensitively on thermal history, since quenching 
and slow cool from high temperature lead to radically different heat capacity and suscepti- 
bility signals 4 . These different approaches offer attractive possibilities to control the physical 
properties in real time. 

The structure of this system creates two different environments for the sodium: Nal which 
lies between adjacent cobalt ions, and Na2 which sits on a lower energy site at the center 
of a cobalt trigonal prism. The delicate balance between these interpenetrating hexagonal 
lattices causes the vacancies on the sodium layer to become attractive at short distances^. 
Long range Coulomb interactions then allow these vacancy clusters to order over long range 
leading to observable superstructure peaks. This Na + ordering buckles the C0O2 layer away 
from the occupied nearest-neighbour sodium sites (Nal) allowing cages to form within which 
the sodium ions are able to rattle. Meanwhile the Coulomb landscape on the cobalt layer 
is modified allowing narrow conduction pathways, and leading to strong correlation of the 
electrons. Both of these properties mean this is a phonon-glass-electron-crystal in agreement 
with studies of Nai.2- x Ca x Co 2 04 and Li .48Na .35CoO^l Dl rt with a 'rattler' site and this 
helps to explain the high thermopower in this system. Variations in the Coulomb landscape 
are expected to lead to the mixed valancy states as observed by NMRPSEl anc i possible 
trapping of spins in the minima. Further observation and understanding of superstructures 
at different concentrations, and of the relevant transitions, will lead to greater insight into 
the overall control exerted on the system by the ionic ordering. 

Whereas Ref. [1] was mainly focused on low temperature ordering, in the range 0.75 < 
x < 0.92, we concentrate here on the evolution of the sodium order as the temperature 
is increased. Unexpectedly we do not observe, in the same composition range, a simple 
disappearance of long-range order at some critical temperature but a rearrangement of the 
sodiums from "square" long range order to striped order at T=285K (as reported iiP, the 
low temperature cell is a slightly deformed square with internal angles close to but not 
equal to 90° that we denote by "square cell" for simplicity). Striped order continues to 
exist, well above room temperature, at least up to T=350K. The reordering of the sodiums 
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is accompanied by a change of the potential exerted by the sodium superstructure on the 
cobalt sites. This has a radical influence on the topology of the conducting paths - which 
have now a unidimensional profile - for the mobile carriers in the Co planes. 

Single crystals of sodium cobaltate were grown using the floating zone technique^ and 
then cleaved to produce high quality samples. Neutron diffraction measurements were car- 
ried out on the SXD diffractometer at ISIS and also on the flat-cone diffractometer E2 at 
HZB. Complementary hard X-ray diffraction on MAGS at HZB was carried out with X- 
ray energy of 12.398keV giving penetration of the sample allowing the bulk ordering to be 
probed. SXD uses the neutron time-of-flight Laue method along with 2tt steradians coverage 
of solid angle by position sensitive detectors to sample large volumes of reciprocal spaced 
This provides an overview of the superstructure peaks across many Brillouin zones which 
were then resolved at higher Q-resolution using E2. Sample conditions on SXD were con- 
trolled using a closed-cycle He refrigerator and by a variable temperature cryostat on E2. 
Three single crystals of nominal composition x = 0.75, 0.92 and 0.78, were investigated at 
different temperatures. The superstructure for x = 0.75 and 0.92 at 150K are reported in 
Ref. m 

Samples of x — 0.75, 0.78 and 0.92 show very similar spectra at T=350K. However, 
at low temperature, the x = 0.78 single crystal showed a striking "hexagon-of-hexagons" 
superstructure (Fig. [TJ, different from that obtained in both x=0.75 and x=0.92 samples^. 
Fig. [I](e) focuses on the hexagon-of-hexagons superstructure pattern that surrounds the 
hexagonal Bragg peak positions on the 1=7 plane for the x=0.78 SXD data. The data show 
a delicate /-dependence (Fig. [l|a, b)) with the superstructure appearing around the Bragg 
positions that agree with h + 2k ^ 3n & Z=even over the range of / measured and then 
around all Bragg positions for odd I values below Z=ll (superstructure peak intensities are 
sometimes weak around the h + 2k ^ 3n positions) and for h + 2k = 3n & 1=11. All 
planes show in-plane modulations of peak intensity towards or away from the (00Z) position 
indicative of buckling of the C0O2 planed 

The superstructure peaks in Fig. [ljfor z=0.78 at T=150K can be indexed on a (^a* x 
Y^b*) grid, and the integrated intensities have been analysed by Reverse Monte Carlo (RMC) 
using simulated annealing in the Canonical Ensemble. The (fictitious) energy E(X) is 
defined in terms of the mean-square deviations of the experimental intensities from the 
calculated structure factors for a given 3D configuration X of all particles in the unit cell. 
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The Metropolis algorithm is used with hopping of a sodium ion from an occupied to a vacant 
site with small displacements of surrounding ions. The temperature is slowly decreased, and 
the configuration at T=0 corresponds to the least-square fit of the integrated peaks. 

A fraction of the peaks can be identified as coming from the square lattice described 
in Ref. DQ and shown in Fig. [2]^a), and these were excluded from the RMC calculations. 
A very robust result from RMC simulations is that the absence of significant intensity 
for most superlattice reflections away from the hexagonal reciprocal lattice points, leads 
naturally to multi-vacancy cluster patterns. For x=0.78 at T=150K we obtain stripes of 
trivacancies condensing into long-range order with cell vectors a' = 5b — a and b' = 5b — 
4a, as shown in Fig. [2|c) with vacancy clusters on the adjacent Na layers being as far 
apart as possible to reduce Coulomb energy giving 3D order. This in-plane stripe structure 
is a simple modification of the square lattice requiring only two shear distortions, of one 
lattice spacing each, to the supercell to reach this end structure (figure [2j.The buckling 
of the C0O2 plane induced by the sodium superstructure has a maximum distortion along 
the c axis of 0.025c; the effect of this buckling on the screening of the Coulomb potential 
is negligible. Fourier transforms of this striped structure are then added to the Fourier 
transforms of the square cell in a ratio of Istr ' pe =1.3. This summation produces remarkable 

-* Square 

agreement with this complex data set (see comparison in fig[j}. Full details of the RMC 
simulations will be reported elsewhere. A number of experiments^^ provide evidence for 
the presence of vacancy clusters, and the formation of these clusters is understood in terms 
of topological constraints and Coulomb repulsiorP^]. This is a first-order effect involving 
energies of a fraction of one eV. Long-range ordering of these clusters is a more subtle 
effect. Refs. Q] and HU agree on the possible occurrence for 0.75< x <0.85 of a number 
of long-range commensurate superstructures with ground-state energies differing by only a 
few meV, which is not significant with respect to second-order effects such as cooperative 
phenomena involving coupling to phonons and mobile carriers in Co layers. We determine 
unambiguously which, among those possible superstructures, the real system chooses. The 
full theoretical interpretation of the results, especially our finite temperature transition, is 
a highly challenging problem. 

The patterns at 350K look quite different from those of Ref. [TJ at T=150 K, with changes 
of symmetry suggesting a phase transition. The superstructure pattern above this transition 
can be labeled on a (|a* x |b*) grid. Figure [3] shows the 1=7, 10, 11 cuts obtained for the 
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£=0.75 samples at T=350K. Again the pattern shows a similar delicate I and in-plane 
dependence. 

X-ray data shows a pronounced hysteresis of over 10K between heating and cooling runs 
(Fig. [4]). Here Q=(0, 1.33,0) is from the square phase in Fig. (2Ta), whereas the peak at 
Q = (0.8, 0, 0) is where the scattering is expected for either of the striped phases in Figs. |2](c) 
and (d). There is a highly unusual peak in the cooling run for Q=(0.8, 0, 0). 

The trivacancy clusters can sit on three possible sites within the stripes in Fig. [2^c). 
Consider a model where the ordering of stripes is long range, but the ordering of trivacancy 
clusters within the stripes is not coherent from one stripe to the next. The translation vectors 
for this phase are a' = 5b + (rj — l)a and b' = 5b + (rj — 4)a where rj — 0, ±1 at random. 
A section of this partially disordered structure is shown in Fig. |2^d). Simulations of the 
scattering intensity from this phase are in excellent agreement with the experimental data 
for x=0.75, 0.78 and 0.92 at elevated temperature, see for example Fig. [3j Therefore, the 
phase transformation is a melting transition from the long-range ordered square or striped 
phases in Figs, pja) and (c) respectively, to the disordered stripe phase in Fig. [2^d). The 
anomalous peak in the intensity of the reflection at Q= (0.8, 0, 0) in the cooling run below 
278K in Fig. gja) can also be explained within this framework. If the disordered stripes in 
Fig. |2](d) transform to the intermediate ordered stripe phase in Fig. |2](c), stronger scattering 
is expected at the same position in reciprocal space, before the reordering to the square phase 
in Fig. |(a) below 274K. 

We note that a different model with quadrivacancy clusters was previously proposed for 
the high-temperature phased However, the model proposed here is a more satisfactory 
explanation due to the simplicity of the phase transformations in fig. [2] and also agrees 
with quantitative RMC analysis. All phases have the same multi-vacancy cluster type 
(trivacancy) and the concentration x in each phase is the same. 

Finally, the entropy associated with the random positions of stripes with respect to one 
another is AS ~ ksNstripes^S, where N str i pes ~ v^A represents the number of stripes (N 
is the number of sites in a Na plane) and there are three different positions within stripes. 
Hence the entropy per particle ^ ~ _strpgs ^ vanishes in the thermodynamic limit. 
Here the disorder is static, but some additional entropy could arise from other static (due to 
defects) or dynamic disorder. Additional disordering of the stripes and longitudinal disor- 
dering of the clusters within stripes may come into play forming something analogous to a 
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smectic phase. The relatively small change in entropy involved in the transition, and any as- 
sociated electronic transitions, are indicated by the lack of any large anomaly accompanying 
the transition. 

This structural transition is accompanied with a topological change of the Coulomb land- 
scape from 2D to ID. Above the transition, there are lines of Co sites at the top of the wells 
(red stripes in Fig. [2]i) which are expected to remain Co 3+ while the conducting holes 
will move along the bottom of the wells (white/yellow stripes). Since the corresponding 
depth, of order 200meV is much larger than the hopping energy t =10meV, we expect at 
room temperature (fcgT ^26 meV) a strong confinement of the carriers along the yellow 
stripes. Patterning of the Coulomb energy could also describe the observation of surface 
superstructure by Pai et al. in which trimers of Na are seen to order in regular arrays^. 
This can be thought as Na sitting on top of the energetically favourable vacancy sites in the 
underlying bulk ordering. Neutron diffraction is not sensitive enough to observe magnetic 
superstructure peaks since this will have an intensity of the order of a hundredth that of 
the magnetic Bragg peak from ref. [TBI The structural transition that we observed through 
neutron diffraction could explain some of the anomalies seen within the 270-290K temper- 
ature range: e.g. broad bump in susceptibility at 285K in x=0.82P', transition in optical 
ellipsometry at 280K in x=0.82 18 , phonon lifetime in infrared conductivity below 295K over 
a range of concentrations^ and anomaly in resistivit y 1 * 20 * 21 !. 

In conclusion, using neutron Laue diffraction and hard X-ray diffraction, we have observed 
square to stripe reordering of the sodium ions near room temperature. As a consequence the 
mobile carriers have to follow restricted paths in the Co planes with a topological change 
from 2D to ID at the transition. This is expected to lead to changes in the topology of the 
Fermi surface, suggesting ARPES measurements to test this prediction. This ID character 
of the conducting path offers an alternative direction in developing new materials with im- 
proved thermoelectric properties. The possibility of producing single domains using electric 
fields or by thin-film deposition should be investigated since the electrical conductivity will 
be higher along the stripes than for multidomain samples, increasing the figure of merit 
for thermoelectric applications. The idea to improve the figure of merit by reducing the 
dimensionality has been suggested a long time agcW Here we find a self-organising system 
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with ID conducting stripes on the nanoscale. 



Electronic address: jonathan.morris@helmholtz-berlin.de 



1 M. Roger, D. J. P. Morris, D. A. Tennant, M. J. Gutmann, J. P. Goff, J. -U. Hoffmann, R. Fey- 
erherm, E. Dudzik, D. Prabhakaran, A. T. Boothroyd, N. Shannon, B. Lake, and P. P. Deen, 
Nature 445, 631 (2007). 

2 J. Braconnier, C. Delmas, C. Fouassier, and P. Hagenmuller, Mat. Res. Bull. 15, 1797 (1980). 

3 F. C. Chou, J. H. Cho, P. A. Lee, E. T. Abel, K. Matan, and Y. S. Lee, Phys. Rev. Lett. 92, 
157004 (2004). 

4 T. F. Schulze, P. S. Hafhger, Ch. Niedermayer, K. Bubenhofer, and B. Batlogg, et al., Phys. 
Rev. Lett. 100, 026407 (2008). 

5 Z. Ren, J. Shen, S. Jiang, X. Chen, C. Feng, Z. Xu, and G. Cao, J. Phys.: Condens. Matter. 
18, L379 (2006). 

6 K. Takahata, Y. Iguchi, D. Tanaka, T. Itoh, and I. Terasaki, Phys. Rev. B 61, 12551 (2000). 

7 I. R. Mukhamedshin, H. Alloul, G. Collin, and N. Blanchard, Phys. Rev. Lett. 94, 247602 
(2005). 

8 H. Alloul, I. R. Mukhamedshin, N. Blanchard, and G. Collin, J. Phys. IV. France 131, 27 
(2005). 

9 M.-H. Julien, C. de Vaulx, H. Mayaffre, C. Berthier, M. Horvatic, V. Simonet, J. Wooldridge, 
G. Balakrishnan, M. R. Lees, D. P. Chen, C. T. Lin, and P. Lejay, Phys. Rev. Lett. 100, 096405 
(2008). 

10 D. Prabhakaran, A. T. Boothroyd, R. Coldea, and N. R. Charnley, J. Crys. Growth 271, 74 
(2004). 

11 D. A. Keen, M. J. Gutmann, and C. C. Wilson, J. Appl. Cryst. 39, 714 (2006). 

12 J. Geek, M. v. Zimmermann, H. Berger, S. V. Borisenko, H. Eschrig, K. Koepernik, M. Knupfer, 
and B. Biichner, Phys. Rev. Lett. 97, 106403 (2006). 

13 F. C. Chou, M. -W. Chu, G. J. Shu, F. T. Huang, W. W. Pai, H. S. Sheu, and P. A. Lee, Phys. 
Rev. Lett. 101, 127404 (2008). 

14 Y. S. Meng, Y. Hinuma, and G. Ceder, J. Chem. Phys. 128, 104708 (2008). 

15 W. W. Pai, S. H. Huang, Y. S. Meng, Y. C. Chao, C. H. Lin, H. L. Liu, and F. C. Chou, Phys. 



Rev. Lett. 100, 206404 (2008). 

S. P. Bayrakci, I. Mirebeau, P. Bourges, Y. Sidis, M. Enderle, J. Mesot, D. P. Chen, C. T. Lin, 
and B. Keimer, Phys. Rev. Lett. 94, 157205 (2005). 

S. P. Bayrakci, C. Bernhard, D. P. Chen, B. Keimer, R. K. Kremer, P. Lemmens, C. T. Lin, 
C. Niedermayer, and J. Strempfer, Phys. Rev. B 69, 100410(R) (2004). 

C. Bernhard, A. V. Boris, N. N. Kovaleva, G. Khaliullin, A. V. Pimenov, L. Yu, D. P. Chen, 
C. T. Lin, and B. Keimer, Phys. Rev. Lett. 93, 167003 (2004). 

S. Lupi, M. Ortolani, L. Baldassarre, P. Calvani, D. Prabhakaran, and A. T. Boothroyd, Phys. 
Rev. B 72, 024550 (2005). 

J. Wooldridge, D. McK. Paul, G. Balakrishnan, and M. R. Lees, J. Phys.: Condens. Matter. 
17, 707 (2005). 

T. Ikeda, and M. Onoda, J. Phys.: Condens. Matter. 18, 8673 (2006). 

L. D. Hicks, T. C. Harman, X. Sun, and M. S. Dresselhaus, Phys. Rev. B. 53, R10493 (1996). 



8 



wavevector Q 1 a* A wavevector Q 1 a* A' 
-4 -2 -4-2 




FIG. 1: (Color) (a, b) Laue diffraction data from SXD for Nao.78Co02 for the Z=ll and 7 planes 
respectively showing the hexagon-of-hexagons around the Bragg peak positions along with the 
corresponding calculations (c, d). The calculation uses phase coexistence of square and stripe 
phases. Detail of the agreement between data for 1=7 (e) and calculation (f). Measured versus 
calculated RMC intensity for the 'square' cell showing the quality of agreement, black line shows 
perfect fit (g). 
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FIG. 2: (Color) Each panel consists of the ionic ordering on adjacent sodium layers (red spheres 
show occupied Nal sites and blue spheres shows Na2 sites as described in the text), in phases 
reported here, separated by the Coulomb landscape on the intervening cobalt layer (white/yellow 
show conduction pathways), a) The square trivacancy cluster structure seen in Nao.TsCoOjH re- 
quires a shear distortion parallel to the a-lattice direction to take it to (b) an unobserved stripe lat- 
tice and then a further shear along b to modify into (c) the stripe structure observed in Nao.7sCo02. 
(d) All three concentrations show the same disordered stripe structure above the melting transition 
at 285K. 
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FIG. 3: (Color) Above 285K, in all three concentrations (x=0.75, 0.78 and 0.92), the long range 
ordered structure (figs. [2^a), (c)) melt into one with disorder between stripes giving the same 
superstructure pattern (fig. [2^d)). (a) 1=11 for Nao.92Co02 at 350K. (b) Z=10 and (c) 1=7 showing 
the I dependence, (d-f) Calculation using the disordered stripe structure. 
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FIG. 4: (Color) Hysteresis behaviour of the ordered phase to disordered stripe transition in the 
intensities of superstructure peaks measured using hard X-ray diffraction on x=0.75. 
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